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Muon-spin rotation spectroscopy (piSR) has been used to 
measure the internal magnetic field distribution in NbSe2 for 
H a i >C H < 0.25 H C 2- The deduced profiles of the supercurrent 
density J s indicate that the vortex-core radius po in the bulk 
decreases sharply with increasing magnetic field. This effect, 
which is attributed to increased vortex-vortex interactions, 
does not agree with the dirty-limit microscopic theory. A 
simple phenomenological equation in which po depends on 
the intervortex spacing is used to model this behaviour. In 
addition, we find for the first time that the in-plane magnetic 
penetration depth X a b increases linearly with H in the vortex 
state of a conventional superconductor. 



Although the bulk of superconductivity research in the 
past decade has been dominated by the study of the high- 
T c cuprates, there remain unresolved issues concerning 
conventional superconductors. In the Meissner state of 
a conventional s-wave superconductor, application of a 
magnetic field H induces pair-breaking and the ensuing 
backflow of quasiparticles gives rise to a nonlinear de- 
pendence of induced supercurrent with respect to the su- 
perfluid velocity [jjj. The solution of the corresponding 
nonlinear London equation gives a penetration depth A 
which changes quadratically with H . This has been con- 
firmed experimentally in Sn, In @ and V 3 Si f|. Theo- 
retical predictions in the vortex state are less developed 
and the ^/-dependence of A in the vortex state of a con- 
ventional superconductor has never been determined ex- 
perimentally. Similarly, the ii-dependence of the vortex- 
core radius po has never been measured in the bulk of a 
conventional superconductor. In a conventional s-wave 
superconductor p ~ £ (Ref. [El), where £ is the coher- 
ence length. Measurements of po bear directly on the 
electronic structure of an isolated vortex and provide a 
necessary comparison for ongoing studies on high-T c ma- 
terials. 

NbSe2 is particularly well suited for a /iSR study of 
the vortex state since the geometry of the vortex lattice 
is well established — thus removing one of the largest ex- 
perimental uncertainties. Scanning tunneling microscopy 



(STM) measurements at the surface and small an- 
gle neutron scattering (SANS) measurements in the bulk 
|8| have produced high quality images of a nearly per- 
fect triangular lattice with long range order. Also the 
coherence length and the magnetic penetration depth in 
NbSe2 are nearly ideal for a /iSR investigation of the vor- 
tex cores. In particular, the large coherence length and 
correspondingly small value of H C 2 (<4 T) imply a large 
signal from the vortex-core regions at moderate magnetic 
fields. 

Most theoretical calculations have modelled the core 
structure using GL theory, which is strictly valid only 
near the superconducting phase boundary. Recently the 
field dependence of the vortex-core radius has been de- 
termined deep in the superconducting state from the mi- 
croscopic theory in the dirty limit, by Golubov and Hart- 
mann @. The vortex -core radius was found to decrease 
monotonically with increasing applied magnetic field due 
to the increased strength of the vortex-vortex interac- 
tions. Although the authors reported good agreement 
with STM measurements || at the surface of NbSe2 at 
T = 0.6 T Cl po was somewhat arbitrarily defined and 
the uncertainty in the measurements was large. The 
results were relatively surprising since NbSe2 is a clean 
superconductor — the ratio of the coherence length to the 
mean free path in the a-b plane is £q/1~0.15 |hJ. In the 
present study, we make a comparison between the vortex- 
core radius measured by /iSR in the bulk of NbSe2 and 
that determined from the spatial dependence of the su- 
percurrent density J s in the dirty-limit microscopic the- 
ory. The results show that the dirty-limit approximation 
is invalid, even at T — 0.6 T c . The if-dependence of po 
is fit to a simple phenomenological equation arising from 
our observation that both X a b and the GL parameter k 
vary linearly with magnetic field. We also report for the 
first time, /iSR measurements of the -ff-dependence of 
X a b in the vortex state of a conventional superconductor. 
Contrary to the H 2 dependence in the Meissner state, 
X a b shows a strong linear 77-dependence in NbSe2 over a 
significant range of applied field. 

In a /iSR experiment the implanted muon samples the 
distribution of local magnetic fields in the bulk of a type- 
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II superconductor in the vortex state. One monitors 
the ensemble averaged muon-spin precession signal. The 
frequency of precession for any given muon is directly 
proportional to the local magnetic field at the muon 
site. Further details regarding the /iSR technique may be 
found elsewhere (see for example Ref. [0). The present 
/iSR study of NbSe2 was performed on the M20 beam 
line at TRIUMF. The single crystal of NbSe 2 used in 
this experiment had a mass of 43 mg and a surface area 
of ^30 mm 2 . The superconducting transition T c and 
H C 2(T = 0) determined from magnetization were 7.0 K 
and 3.5 T, respectively. For details regarding the sam- 
ple growth see Ref. Q. The sample was mounted with 
its c-axis parallel to the applied field and beam direction. 
/iSR spectra were recorded under conditions of field cool- 
ing in a He gas flow cryostat. A cup-shaped veto counter 
was used to suppress unwanted background signal from 
muons which missed the sample |l3| . 

In Fig. 1 the Fourier transforms of the muon preces- 
sion signal in NbSe2 are shown for different fields at 
T = 0.33 T c . The real amplitude of the Fourier trans- 
forms is a good representation of the internal magnetic 
field distribution from the vortex lattice convoluted with 
small nuclear dipolar fields. Note the small peak near 
zero in the top panel of Fig. 1. This is due to a small 
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FIG. 1. The Fourier transforms of the muon spin precession 
signal in NbSe2 after field cooling to T = 0.33 T c in magnetic 
fields of H — 0.1, 0.35 and 0.7 T. The average magnetic field 
of the residual background signal is denoted as Btkgd- 



(2%) residual background signal due to muons which miss 
the sample. Lineshapes have been renormalized to the 
same maximum amplitude. The sharp features expected 
for a perfect triangular vortex lattice, such as the Van 
Hove singularity at the saddle point, are obscured partly 
by broadening effects of the finite Fourier transform and 
by flux-line lattice disorder. Nevertheless there is a clear 
high-field cutoff observed in the /iSR lineshape originat- 
ing from the finite size of the vortex cores. The effect 
of increasing H on the high-field cutoff is clearly seen in 
Fig. 1. At all of the magnetic fields the signal-to-noise 
ratio of the high-field tail is so large that one can unam- 
biguously extract the vortex-core radius. 

In order to test the strength of the pinning forces on 
the vortex lattice the sample was cooled in an applied 
field of 0.5 T to 2.3 K after which the field was decreased 
by 7.5 mT. In our previous studies of YBa2Cu30y_5 @] 
we found under similar conditions that the background 
signal shifted to a lower frequency in response to the 
change in the applied field, while the signal corresponding 
to muons stopping in the sample remained unchanged. 
In NbSe2 both signals shifted — indicating that the vor- 
tex lattice is considerably more weakly pinned than in 
YBa 2 Cu 3 7 _ 5 . 

The /iSR spectra were fit in the time domain where 
there are no complications associated with fitting finite 
Fourier transforms [ jL5| . The distribution of muon pre- 
cession frequencies from the vortex lattice was modelled 
with a theoretical field distribution generated from a GL 
model jl6). The local field at any point in the a-b plane 
is given in a suitable approximation by |17j], 
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Here Bq is the average magnetic field, G are the recip- 
rocal lattice vectors, b = B /B C2 , £ a {, is the GL coherence 
length and K\ (u) is a modified Bessel function. The cut- 
off factor u K\ (u) accounts for the finite size of the vortex 
core — whereas in the London model B(p) di verg es loga- 
rithmically as p^0. Recently, Yaouanc et al. |L7[ showed 
that u K\ (u) is a good approximation of the cutoff factor 
determined from the exact numerical solutions of the GL 
equations |l8) at low reduced fields b. 

The theoretical muon polarization function was gen- 
erated by assuming the field profile of Eq. (Q) and then 
multiplying by a Gaussian relaxation function e~ a 1 / 2 
to take into account any residual disorder in the flux- 
line lattice and the contribution of the nuclear dipolar 
moments to the internal field distribution. The residual 
background signal was fit assuming a Gaussian broad- 
ened distribution of fields. All fitted parameters were 
treated as independent variables. From the fitted values 
of cr, the RMS deviation of the vortices from their ideal 
positions in the triangular lattice was determined to be 
less than 3% of the intervortex spacing over the entire 
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FIG. 2. The magnetic field dependence of (a) X a b(H) and 
(b) k'(H) = X ab (H) / p (H) in the vortex state of NbSe 2 at 
T = 0.33 T c (open circles) and T = 0.6 T c (solid circles). The 
solid line fits are described in the text. 

field range studied. This small disorder is consistent with 
STM and SANS imaging experiments on NbSe2. 

The magnetic field dependence of X a b is shown in 
Fig. 2(a). Contrary to the Meissner state, a linear- H 
dependence is observed in the field range studied. A 
fit to the linear relation X a b(H) — A Q i,(0)[l +{3h], where 
h = H/H c2 (T), gives A ab (0) = 1323 A and (3 = 1.61 
at T = 0.33 T c and A a6 (0) = 1436 A and (3 = 1.56 at 
T = 0.6 T c . We note that d[AX/X(0)}/d{H/H c2 ) is con- 
siderably weaker than for YBa 2 Cu3 06.95 (Ref. JT^j), in 
which there is strong evidence for line nodes in the su- 
perconducting energy gap function. 

We define an effective vortex-core radius po to be the 
distance from the vortex center for which the supercur- 
rent density J s (p) reaches its maximum value. J s (p) was 
obtained from fits of the data to Eq. (|l|) and the Maxwell 
relation J(p) = VxB(p). In Fig. 3 /xSR measurements 
of po are shown as a function of H at T = 0.33 (open 
circles) and 0.6 T c (solid circles), along with the STM 
measurements (open squares) of Ref. . The smaller er- 
ror bars and reduced scatter in the pSR data reflects the 
statistical improvement of a /iSR experiment which sam- 
ples a large number of vortices in the bulk of the crystal, 
as opposed to STM which averages the radius of a few 
vortices at the surface. The dashed curve drawn through 
the STM results comes from tunneling current I(p) pro- 
files calculated from the Usadel equations, as explained 



in Ref. |)| . To generate J s (p) profiles from Usadel's dirty- 
limit theory we extend the work of Ref. |j| to include the 
self-consistency equation for the vector potential A(p). 
In cylindrical coordinates the equation of motion is [pL9|| , 

- (p^-) = fi~ 2 A 2 sin 9 cos 9 - A cos 9 + lu sin 9, (2) 
pap \ dp J 

where 9 parametrizes Usadel's normal (G = cos 9) and 
anomalous (F = sin9) Green's functions |20[ ], A(p) is the 
order parameter, u> — (T/T c )(2l + 1) is the Matsubara fre- 
quency and k = [47t 5 /7('(3)] 1 / 2 k (where k = A/£). Equa- 
tion (g) is supplemented by the following self-consistency 
equations, 

A In (T/T c ) = -2(T/T C ) £ [A/w - sin 9} (3) 

UJ 

Js{p) = ~ (j-pPA) = 16^- 2 (T/T c ) A ^ sin 2 0, (4) 

and the boundary conditions for singly quantized vortices 
with a Wigner-Seitz cell radius p s = (Qq/ttH) 1 / 2 , 

A(0) = 0(u, 0) = 0, A'(p s ) = 9'(uj, p s ) = (5) 
A(p^0)^~R/p, A( Ps ) = Q. (6) 

Equation (^|) subject to these boundary conditions was 
solved numerically for 9(to,p) starting with the ini- 
tial trial potentials A(p) = Aq tanh(p) and A(p) — 
R (l/p— p/ Ps). Improved values of A(p) and A(p) were 
obtained by including the self-consistent conditions (|^) 
and (|]). The parameter k in equations (||) and (||) was 
determined from fits to the data. We found k to be nearly 
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FIG. 3. The magnetic field dependence of the vortex-core 
radius in NbSe 2 determined by STM [| at T = 0.6 T c (open 
squares) and by ^tSR at T = 0.33 T c (open circles) and 
0.6 T c (solid circles). The solid curves are from calculations 
of supercurrent density J s (p) profiles using H C 2(0) = 3.5 T, 
ff c2 (0.33 T c ) = 2.9 T and H c2 (0.6 T c ) = 1.9 T from magne- 
tization. The dotted curves through the /iSR data are from 
Eq. (§. 
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temperature independent at all fields studied, which is 
consistent with the original definition of k near T c in 
GL theory The deduced values of po were not very 
sensitive to n. In Fig. 2(b) we define k' — \ a b/po, where 
k'{H) = \Mk{H) - 1.98 and k'(H) = 1.12k(H) - 1.26 
at T = 0.33 T c and T = 0.6 T c , respectively. At both 
temperatures k! (and hence k) increases linearly with H. 
Fitting to the linear relation n'{H) = k'(0)[1 +jh], we 
obtain k'(0) = 6.9 and 7 = 9.5 at T = 0.33 T c and k'(0) = 
5.1 and 7=10.2 at T = 0.6 T c . 

Figure 4. shows the theoretical J s (p) and A(p) profiles 
together with the J s (p) profile obtained from experiment 
for a particular T and H. The vortex-core radius taken 
from the J s (p) profiles of the dirty-limit theory are shown 
as solid curves in Fig. 3. Not surprisingly there is poor 
agreement with the fiSK data at T = 0.33 T c , where ther- 
mal smearing of the bound states in the vortex core is 
negligible. Contrary to the STM results ||, however, 
there is also poor agreement at T = 0.6 T c — suggesting 
that the dirty-limit theory does not adequately describe 
the shrinking of the vortex-core radius with increasing 
H. In the STM experiment it was necessary to arbitrar- 
ily define po as the radius in which the tunneling current 
decreased to 36% of 7 max - I m i n and A(p)/A( j o s ) = l/v / 2- 
However, we find that the J s (p) profiles generated from 
the dirty-limit theory do not peak exactly at a radius 
corresponding to A(p)/A(p s ) = l/\/2 for all values of 
T/T c and H/H C 2- Thus our definition of po should pro- 
vide a better description of the true -ff-dependence of the 
vortex-core radius. The STM data may also be influenced 
somewhat by the discontinuity in the energy spectrum of 
the vortex cores which occurs at the sample surface. It 
has been suggested that the effect may be an enlargement 

of P0 m. 
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FIG. 4. The current density and order parameter for 
T = 0.6 T c and H/H c2 (0.6 T c ) = 0.053. The solid curves 
are from the dirty-limit microscopic theory while the dashed 
curve is the supercurrent density from a jj.SK measurement of 
the field distribution in NbSe2. 

On the other hand, the /iSR results fit well (see dotted 
curves in Fig. 3) to the phcnomcnological equation, 



where po(0) = A a b(0)/«; / (0). From the fits to X a b(H) and 
k'(H), p (0) = 191 and 282 A at T = 0.33 and 0.6 T c , 
respectively. For a triangular vortex lattice H = Bo = 
3<&o/\/2£ 2 , where L is the intervortex spacing. Thus 
for a given temperature Eq. (Q) may be rewritten as a 
function of the distance between vortices. 

In conclusion, we have determined that the magnetic 
penetration depth A a f, shows a linear magnetic field de- 
pendence in the vortex state of NbSe2 . The linear term is 
considerably weaker than that determined previously for 
YBa2Cu30e.95. Also, we find that the vortex-core radius 
po which shrinks with increasing field is not adequately 
described by the dirty-limit microscopic theory, but does 
obey a simple phenomenological equation involving the 
intervortex spacing. Our results imply that the conven- 
tional GL equations with field- independent length scales, 
are not applicable deep in the superconducting state. 
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